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Abstract. The news function providing some relevant information about angular distri- 
bution of gravitational radiation in axisymmetric black hole collisions at the speed of light 
had been evaluated in the literature by perturbation methods, after inverting second-order 
hyperbolic operators with variable coefficients in two independent variables. More recent 
work has related the appropriate Green function to the Riemann function for such a class 
of hyperbolic operators in two variables. The present paper obtains an improvement in 
the evaluation of the coefficients occurring in the second-order equation obeyed by the 
Riemann function, which might prove useful for numerical purposes. Eventually, we find 
under which conditions the original Green-function calculation reduces to finding solu- 
tions of an inhomogeneous second-order ordinary differential equation with a non-regular 
singular point. 
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The non-linear Einstein equations [1] ruling the gravitational field are so complicated 
that no exact solution has been obtained so far without making a number of simplifying 
assumptions. On the other hand, in the physically more relevant case of isolated gravitating 
systems, which are time-dependent, no simplifying assumption can be made apart from 
axisymmetry. Within this framework, many efforts have been devoted to analytic and 
numerical investigations of gravitational radiation produced in axisymmetric black hole 
collisions at the speed of light [2-5] , since such events (although unlikely) are expected to 
lead to the largest amount of gravitational radiation ever studied at theoretical level. 

The angular distribution of gravitational radiation is described in part by the news 
function [6], and the work by D'Eath and Payne [2-5] has shown that, in the above events, 
such a function can be obtained at second order in perturbation theory provided that one 
is able to solve a set of inhomogeneous hyperbolic equations taking eventually the form 
(m and n being integers) 

^m,nX{Q,r) ^ H{q,r). (1) 

With the notation in Refs. [2-5], q and r are the independent variables, the function H is 

related to the source term in the original set of equations, and the operator Crn,n reads 

= -(2^/2 + 32,)^ + 4,^^ + 64|1 + 4(n + 

d 

- 16n— + n'^ - . (2) 
or 

The inverse of the operator Crn,n is an integral operator whose kernel is equal to the Green 
function Gm,niQ,r; qo,ro) of Eq. (1). 

The work in Ref. [7] has however pointed out that, after reduction of Eq. (1) to 
canonical form through the introduction of suitable new variables: 



t = \J 1 + 16gV2, (3) 



x = r + \og[ ^ ) - - 4, (4) 
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,.. + log(^)+^-4, (5) 

Eq. (1) can be solved for x with the help of a standard integral formula [8] which involves 
the Riemann function. This is a valuable tool in the theory of hyperbolic equations in two 
variables with variable coefficients, but unfortunately it has not been much exploited (to 
our knowledge) in the literature on gravitational physics. Following Ref. [7] , we consider 
the operator 

having defined w as the function oi x — y such that 

w^-l x-y 

we 2^ =68, (7) 

from which 

t=i±^. (8) 

1 — W 

r / / 1 (1 - + l)2(2n + 1 + t) 
a{t) = a[t{w{x - y))] = (^4 + 4^2 _ 1) ^' (9) 

,/N ,r/ / 1 (t + l)(t-l)^(2n + l-t) 

c(t)=c[tM.-y))]='"'^-"''';-')''!+^;. (11) 

w L V V 256 {t^ + 4*2-1) ^ ^ 

The coefficient of the highest-order derivative in Eq. (6) is constant because 

^m,nX = fLx = H, 
where / is a function evaluated in Ref. [7] , and hence we study the equation 
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eventually, whenever / does not vanish. The (formal) adjoint of the operator L is then 
equal to 



where 



A = -a, 
B = -6, 



^ _ da db 
dx dy 



(12) 

(13) 
(14) 

(15) 



Let now (^,?7) be coordinates of a point P such that characteristics through it intersect 
a curve 7 at points A and -B, AP being a segment with constant y, and BP being a 
segment with constant x. As a function of x and y, the Riemann function R satisfies the 
homogeneous equation 

L^R = Q, (16) 

the boundary conditions 



dR ,^ dR 

— = bR on AP, — = aR on BP, 

ox oy 



(17) 



and is equal to 1 at P. It is then possible to express the solution x oi Eq. (1) in the form 
[7] 



x{P) = l[x{A)R{A) + x{B)R{B)]+ f 
^ Jab 



Rdx A „ IdR, 

2a^ + i^^-2a^>^ 



dx 



Rdx ( ^ ldR\ 

2a? + r«-2a^)''. 



dy 



+ 



/ / R{x,y\i, 
J Jo. 



v)- 



-H{x,y) (t-l)2(t + l)2 
256 (t^ + 4^2 _ 1) 



{x, y)dx dy, 



(18) 



Vt being a suitable bounded domain. 

The main technical difficulty in Eq. (16) results from the derivatives occurring in 
the coefficient C; since w{x — y) is found only implicitly from Eq. (7), and hence t, a, b, c 



from Eqs. (8)-(ll), it would be helpful not having to take numerical derivatives in (15) 
of a function which is only found numerically itself. The present note solves this specific 
problem by pointing out that, in Eq. (15), one has 



da da dt dw 
dx dt dw dx ' 



(19) 



where, from Eq. (9), 



da 
di 



1 
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(4^3 + 6(n + l)t^ + Ant - 2{n + 1)) 
(t4 + 4t2 - 1) 



4t(l -t){t + 1)2(^2 + 2)(2n + 1 + t) 
(^4 + 4^2 - 1)2 



while Eq. (8) yields 



dt 



dw (1 — w)"^ 

Moreover, from the logarithm of both sides of Eq. (7) we find 



log('u;) + 



w 1 (x — y) 
2 ~ 2^ ~ 8 



The derivative with respect to x of both sides of Eq. (22) yields eventually 



dw _1/ 1 
dx 4 V w 



-2 



(20) 



(21) 



(22) 



(23) 



Equations (19)-(23) yield a complete algorithm for the evaluation of the first derivative in 
C, i.e. 

da 1 da w"^ , , 

di " 2 at (1-^2)2' ^ > 

where Eqs. (7), (8) and (20) should be inserted for the purpose of numerical analysis. An 
entirely analogous procedure holds for |^ in Eq. (15), i.e. 



dh db dt dw 
dy dt dw dy ' 



(25) 
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bearing in mind Eq. (10), and that 



which leads to 



with 



dw 
dy 



db 1 db 



dw 
dx ' 



w 



db 
di 



1 

16 



dy 2dt{l-wY' 



(4^3 - 6(n + 1)^2 + 4nt + 2{n + 1)) 
(t4 + 4*2 - 1) 



4:t{t + l){t-iy{t^ +2){2n+l-t) 
(^4 + 4*2 - 1)2 



(26) 



(27) 



(28) 



Eventually, the coefficients A, B and C in the operator are all expressed as functions 
of w = w{x — y) with the help of previous equations as follows: 



A = 



w(l + n — nw) 



4(1 + 2w- 2^2 + 2^3 + ^4) 



B 



w'^{n{—l + w) + w) 
4:{1 + 2w - 2w^ + 2w^ + w^) 



(29) 



(30) 



64(1 + 3w + 3«;^ + w^fC = (m^ - n^){l + 4t(; + Bw^ + Sto* + 4t/;^ + w'') 



2 „2> 



.5 , „..6> 



+ 4(n + 1)(1 - - -w^ + w^) 



(31) 



where we have used Eqs. (8)-(ll), (13)-(15), (20), (24), (27) and (28). We have solved 
numerically Eq. (7) after re-expressing it in the form 



8 log('u;) -\- A[w ) — X — y, 



w 



(32) 



for a; — y in the closed interval [0.8, 10], at steps of 0.1. The variable w is then monotonically 
increasing and ranges from 1.05 through 1.83. The corresponding values of A, B and C 
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have been obtained upon insertion of w into Eqs. (29)-(31). For example, at the initial 
point X — y = 0.8, the coefficients A, B and C read 



A = 0.059 - O.OOaOn, (33) 

S = 0.065 + 0.0031n, (34) 
C = -0.0077 + 0.0038(m2 - n^) - 0.0077rz, (35) 

respectively. 

Work is now in progress in applying such formulae to the numerical evaluation of 
the Riemann function itself [9]. Hopefully, this will lead to more powerful tools for the 
investigation of gravitational radiation in the few cases where it is expected to be very 
rich, i.e. axisymmetric black-hole collisions at the speed of light [2-5] . 

Meanwhile, we are also making progress on alternative approaches to finding solutions 
of the equation Lx = y = H that we started with. To exploit both Fourier transform 
theory and the fact that the coefficients in the operator L depend on t and hence, eventually, 
only on the diflFerence x — y,we give up the canonical form of the operator L by introducing 
the independent variables 

x.^. y.^. (36) 

Our inhomogeneous hyperbolic equation is then turned into 

Lx{X, Y) = 4H{X, Y) = h{X, Y), (37) 

where 

dX^ dY^ ^ 'dX ^ 'dY ' ^ ^ 

with coefficients depending only on y. It is therefore more convenient to consider the 
Fourier transform of x with respect to X (assuming that it exists), i.e. 



1 

x{k,Y)^-= / x{X,Y)e-''^''dX, 



(39) 
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with inversion formula 



1 1"°° 

V^TT J-oo 



(40) 



If also the right-hand side of Eq. (37) admits Fourier transform h{k, Y) with respect to 
X, we therefore obtain the second-order equation 



^ + 2(5 - A) A + _ 4C - 2ik{A + B)) 



By virtue of Eqs. (9), (10), (13) and (14) we find 



x{k,Y) = -h{KY). (41) 



B-A 



1 (n + 2)t(l -t^; 
8 (^4 + 4^2-1) 



l (t^-l)(t^ + 2rz+l) 
" 8 (t4 + 4t2-l) ^ 

and we re-express Eq. (32) in the form 



(42) 



(43) 



41og(Tij) + 2 [w--^] = Y. 



(44) 



Note now that, for each fixed value of k, Eq. (41) may be viewed as an inhomogeneous 
second-order ordinary differential equation for x- On denoting by M the operator in square 
brackets on the left-hand side of Eq. (41), and by xi ^^'^ X2 two linearly independent 
integrals of the homogeneous equation Mx = 0, one can use the method of variation of 
parameters, according to which the general solution of Eq. (41) reads 



X = XiXiiY) + A2X2(1^) + viiY)xiiY) + V2iY)x2iY), 



(45) 



with Ai and A2 constants, while vi and V2 are chosen in such a way that 



^1X1 + V2X2 = 0, 



(46) 



vWi + vW2 = -h = jiY), 



(47) 



the prime denoting derivative with respect to Y. The system expressed by Eqs. (46) and 
(47) is solved by 

X2{Y)j{Y) 



vi{Y) = -l 



V2iY) 



W(Y) 

Xi{Y)j{Y) 
W{Y) 



-dY, 



dY, 



(48) 
(49) 



where W is the standard notation for the Wronskian of xi ^ind X2, i-e. 

W{Y) ^ XiiY)x2iY) - X2{Y)x[{Y). (50) 

In the course of investigating the singular points of Eq. (41), the technical difficulty results 
from the fact that its coefficients are ratios of polynomials in the t or w variable, but not 
in the Y variable itself by virtue of Eq. (44). A way out is obtained by turning Eq. (41) 
into an ordinary diflFerential equation with w taken as the independent variable. For this 
purpose, we exploit Eq. (44) to find 

-2 



dw 



1 



dY 2\^w 



1 



and hence 



^+^l(-)^+P2(«;) 



X{k,w) = a{w), 



(51) 



(52) 



where 



giiw) = 



1 + 



w 



+ 4:{B-A){l + 



1 



w 



(n + 2) 



(l + «;)2(l-w2) 



w{l + w) 



(1 + (1 + ^/^)w + w;2)(l + (1 - ^/l)w + w"^) ' 



w 



(53) 



w 



^2(w)=4 1 + - {k^ - AC -2ik{A + B)) 



(1 + ^1^2 1 



IGu'^ 



(m^ — m?)- 



1 + (1 + Vh)w + w2)(l + (1 - V^)w + 



(1 - w2 - -w^ + w^){l + wf 
(1 + (1 + ^/h)w + w'^Y{l + (1 - \/l)w + w'^Y 

((n + 1)(1 + w^) - 2nw){l + 

2 w^{l + (1 + \/l)w + t(;2)(l + (1 - s/l)w + w"^) 
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+4(rz+l) 
ik 



(54) 



1^^ 



w 



a{w)=4ll + - j{Y{w)). (55) 



Equation (52) has therefore a non-regular singular point at w = 0, in that the coefficient 
g2{w) has a fourth-order pole therein. This corresponds to the point x = — oo in the 
original variable defined in Eq. (4). Moreover, the polynomial p{w) = 1 -|- (1 + y/5)w + viP 
has real roots equal to 



w^,2 ^ + V5) ± ^^(l + \/5). (56) 

At these points, g\{w) has a first-order pole, and g2{w) has a second-order pole. The last 
pole occurs for gi{w) at = —1, and it is of first order. No further poles occur at finite 
real values of w, since the polynomial q{w) = 1 -|- (1 — \/5)w + w'^ has the complex conjugate 
roots 



By virtue of (54), Eq. (52) does not possess normal integrals. These are meant to 
be integrals admitting the factorization e^u [10], where is a polynomial in ^ and u 
solves an equation with a Fuchsian singularity at the origin. The point w = would be 
therefore an essential singularity through the occurrence of the factor e^. In our case, on 
re-expressing gi and g2 in the form 

w[l + w) w 

g2{w)=4 ^ — k + — + — , 

with obvious meaning of the functions Fi,F2,F3 by comparison with (53) and (54), the 
ansatz for normal integrals yields the following second-order equation for u: 



+(2n' + ^i)A + (n" + n'VpiO' + ^2) 



dw'^ dw 



u = e-^a. (57) 
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Since u should obey an equation with regular singular point at the origin, we try to choose 
in such a way that the coefficients of w~'^ and vanish in Eq. (56). The former task 
is accomplished by choosing 

n' = 2ik{i + - 

\ w 
which yields 

n" + —(fi + 7-^-^ ) + ^ = A (^3 + 2i/c(l + wfFi + 4ikw) + 0{w-^), (58) 
w \ [1 + w) J w-^ ^ ' 

where the term in round brackets multiplying does not vanish. This is why normal 
integrals cannot be found in our problem. 

We have therefore to look for x the most general form [10] 

oo 

X{k,w) = w'^ ^ aiw\ It; e]0,i(;i[, (59) 

l= — oo 

where 7 accounts for the multi-valuedness of the solution (any integer part being absorbed 
in the Laurent expansion valid in the open interval w e]0,wi[). This expansion is now 
inserted into the homogeneous equation associated to Eq. (52), i.e. 



+ ^l(w)— +fi'2(w) 



dxir dw 



X{k,w)=0, (60) 



since the knowledge of two linearly independent integrals of Eq. (60) is sufficient to solve 
Eq. (52) by exploiting the method of variation of parameters previously described. For 
our purposes we consider Laurent expansions of gi and g2 in the form 

00 h ( \ 

g,{w)= J2 aV = ^^, we]0,w,[, (61) 

s=— 00 



92{w)= b^w^ = ^, we]0,w,[. (62) 
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We therefore get from Eq. (60) the condition 



p= — oo 



^ j{r,i,p,k) 



ai 



.l=—oo 



njP = 0, 



(63) 



having defined 



J{r, l,p, fc) = (7 + 0(7 + ^ - l)Si,p+2 + (7 + + bUik)- (64) 

Equation (63) leads to the infinite system of equations 



XI JiT,l,P,k)ai = 0, 

l= — oo 



(65) 



for all p = — oo, ...,+oo, where the coefficients and occurring in J are all known. 
Indeed one finds (cf. (53) and (54)) 



l-^-2w-nw^ + 0{w^) 



w 



(66) 



92{w) = 



(-1 + 4A;2 - n) + 2(1 + 8k'^ + n - ik{l + n))w + 0(«;^ 



if w e]0,t(;i[. By virtue of (64), Eq. (65) may be cast in the form 



1— — 00 



having defined 



(67) 



(68) 



G^ip,l) = 



1 



(7 + p)(7 + p- 1) 



(69) 



where the general form of the coefficients and is found to be, in agreement with (53) 
and (54), 

a'l =nai^i + a2,i, (70) 
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6° = m^as^i + n^a^^i + na^^i + cce,/. (71) 

The numerical coefficients ak,ii for all /c = 1, 6, can be inferred from (66) and (67). For 
example, one finds 



q;i,i = --, a2,i = "1,2 = -2, 0:2,2 = -2, cki.s = 7, 0:2,3 = 12. 

The polydromy parameter 7 is found, at least in principle, by requiring that the homo- 
geneous linear system expressed by Eq. (68) should have non-trivial solutions. Such a 
problem is currently under investigation. Hopefully, an intriguing link between black hole 
physics and the theory of infinite determinants [10] will be found to emerge. 
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